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A novei statisticai modei for the cooperative binding of monomeric iigands to a linear iattice is deveioped to 
study the interaction of ionic surfactant moiecuies with flexibie poiyion chain in diiute soiution. Eiectrostatic 
binding of a ligand to a site on the polyion and hydrophobic associations between the neighboring bound 
ligands are assumed to be stochastic processes. Ligand association separated by several lattice points within 
defined width is introduced for the flexible polyion. Model calculations by the Monte Carlo method are carried 
out to investigate the binding behavior. The hypothesis on the ligand association and its width on the chain are 
of importance in determining criticai aggregation concentration and binding isotherm. The resuits are reason¬ 
able for the interpretations of several surfactant-flexible polyion binding experiments. The implications of the 
approach are presented and discussed. 
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1. Introduction 

Complex formation of surfactant molecules with polymer chain is one of the most important and 
attractive subjects in coUoid and polymer science. It is also useful in various fields of applied chemistry, 
such as pharmaceutical chemistry, food science, cosmetic manufacture, and so on. Extensive studies have 
been carried out on interaction of ionic surfactants with oppositely charged linear polyion, in particu¬ 
lar QHs). 

One of the interesting subjects to be revealed in this field is the cooperative nature of surfactant 
binding due to the hydrophobic interactions among bound surfactant molecules on the polyion chain. 
To know about it, it is essential to understand the behavior of binding isotherms as weU as the critical 
aggregation concentration (CAC) which characterizes the onset of surfactant binding to polyion (4|[3. The 
cooperative nature is due to the side-by-side hydrophobic interactions (association) of the aliphatic tail of 
the surfactant molecules bound to the polyion chain. 

Concerning the theoretical approaches, simple cooperative site binding theory based on the Ising 
model has been widely used to successfully analyze the binding isotherms in most cases of dilute solu¬ 
tions of stiff polyions (6][7l. In our previous works, the effects of steric hindrance and of intrapolymer 
cooperative interaction between bound hgands across an unoccupied binding site, named skip-binding, 
are investigated by expanding the simple Ising model. The effect of steric hindrance due to head group 
size should be taken into account in some systems (8]|9l. Data fittings using matrix method to this modified 
binding model show that this approach is available to the analysis of binding isotherms to stiff polyions 
(ini. Furthermore, it is important to consider the effect of polyion flexibility on the surfactant binding to 
understand several experimental studies I11H13I . 
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However, there seem to be still at least two shortcomiitgs in these models at the moment. Firstly, the 
aliphatic tail of the bound surfactant molecules has some degree of freedom with its binding orientation: 
this should lead the neighboring interactions between the bound surfactants to be stochastic process 
essentially. Secondly, the hydrophobic association between the bound surfactants separated by several 
sites could substantially take place in the ionic surfactant-flexible polyion system. These two new points 
of view should be taken into account in the site binding model treatment. 

In the improved model, the polyion molecule, on which bound monomeric surfactant molecules as¬ 
sociate with each other, is also assumed to be a linear lattice chain. Here, electrostatic ligand binding 
to a lattice point and the hydrophobic association between bound ligands are assumed to be stochas¬ 
tic processes according to the free energy change. In addition, the associations between bound ligands 
are acceptable, when their separation is within the given range on the linear lattice. This expansion of 
the binding model is necessary for a sufficiently flexible polyion chain. The calculation procedure using 
random numbers is employed to obtain the binding isotherms of the system. 

In this paper, the typical calculations of the parameter dependence are attempted for the interpreta¬ 
tion of the binding isotherm of ionic surfactants to the flexible polyion with standard Monte Carlo (MC) 
algorithm. The results of this simple one-dimenslonal model calculation show an apparent increase of 
the surfactant binding, and is reasonable for the interpretation of the characteristics observed in the ex¬ 
perimental studies. The expansion of the scheme on the association of a pair of the hgands is essential 
for the binding to the flexible polyion. The present approach for a new lattice model should be feasible 
to analyze the experimental results of the system of ionic surfactant-flexible polyion. 


2. Theory 

2.1. The model and parameters 

In this model, a ligand is not necessary to interact with the ligands on the adjacent binding sites, 
differently from the ordinary Ising models Including our previous ones. It is based on the assumption 
that the electrostatic interaction permits the ligand bindings in spatially different orientations around 
the polyion, although the model does not explicitly include the electrostatic interaction. Another signifi¬ 
cant point of the model is the introduction of the association across the multiple sites. The characteristic 
parameter association width n, mentioned later on. Is introduced to the model to reflect the persistence 
length of a polyion. Namely, the more flexible fhe polyion backbone is, the larger n the polyion has. 
These assumptions enable the association between ligands, separated by several (occupied and/or non- 
occupied) binding sites on the polyion. This brings a partially two- or three-dimensional character into 
the one-dimensional lattice, and improves the behavior of ligand associations which is more suitable for 
analysis. 

A linear chain consisting of n equivalent charged lattice points (binding sites) is regarded as a model 
of a flexible polylon. A monomeric ligand (oppositely charged ionic surfactant molecule) can bind to one 
of the sites with the binding free energy change Acb, where the intrinsic binding free energy change is 
Ae°. The main part of the change Ae° is due to electrostatic energy change between a charged site and 
an ionic head group of the surfactant. Bound ligand is capable of stochastically associating with other 
bound ligands with the association free energy change Ae^ primarily due to the hydrophobic interac¬ 
tions between the aliphatic tails. This association is assumed to cover the bound ligand on the n-\h neigh¬ 
boring lattice points on both sides due to polyion chain collapse. The free energy change contains the 
contribution from the polyion deformation. The number n, dependent on polymer configuration, poly¬ 
mer flexibility, and surfactant size, is called ‘association width’ in this manuscript. In the case ol n^2, 
the association is possible between the bound ligand separated by lattice point(s). One ligand can simply 
associate with two ligands on both sides of the polyion chain in the present model, as shown in figurej^ 

The pair association between the bound ligands is chosen to be the nearest neighbor within the asso¬ 
ciation width n. In this model, electrostatic hgand binding without hydrophobic association is permitted 
even within the association width. This association does not prevent the associations by other ligands 
between the pair (see figurej^. In our previous work, cooperative ligand association across only one un¬ 
occupied site (skip-binding) was introduced for the surfactant association on the relatively stiff chain ITOl . 
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Figure 1. Schematic representation of the surfactant binding to the sites on the flexible polyion in the 
present model. The possible pair associations between bound hgands in the case of ‘association width’ 
n = 4 are indicated by thick arrows between aliphatic tails of surfactant molecules. 

The assumption in the present model is expanded more suitably for the flexible polyion chain, although 
the formation of very long loops is ignored and the change of the model makes the analytical calculation 
difficult. 

The other quantities and the parameters used in this manuscript are defined as, Cf, free ligand concen¬ 
tration, 6, binding fraction of sites (degree of binding), K, intrinsic binding constant, and u, cooperativity 
parameter where a -llu was used as a cooperativity parameter in our previous papers. 

In the present model, five (one ligand-free and four bound) states of a binding site are assumed on 
each lattice point. Their symbols and definitions are 


f, 

bo, 

be, 

bR, 

b2, 


there is no bound ligand, i.e., the site is ligand-free, 
there is a bound ligand which does not associate with other ligands, 
there is a bound ligand which associates with a ligand on left side, 
there is a bound ligand which associates with a ligand on right side, and 
there is a bound ligand which associates with two ligands on both sides. 


A bound ligand without associated ligands within the width n is in state bo. If the ligand has associable 
ligands, its state is stochastically determined according to their energy change in the MC trials. The states 
and their transitions are schematically illustrated in figurej^ 

The fraction of sum of bound states (binding fraction) 6 is represented as. 


0 = [bol + [bL] + [bR] + [bzl 


( 2 . 1 ) 


where the symbol with brackets means their fraction on the polyion and obviously [bR] = [bR] since 
the binding with left-hand and right-hand ligands is equivalent. The relations of K and u with these 
concentrations and the above mentioned free energy changes are expressed as. 


■|^ = exp(-A£-°/fcBr), 


( 2 . 2 ) 



(2.3) 


a 


respectively, where is the Boltzmann constant and T is absolute temperature. 
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Figure 2. Schematic representation of five states (f, bo, bL, b^ and b 2 ) of a binding site (small closed circle) 
on the polylon chain with surfactant and their associations (small arrows). Large arrows with a fraction 
indicate the directions of the state transitions and their probabilities of the MC trial, respectively. 


The formation of nano- and microparticles is one of the interesting subjects of the surfactant-polyion 
(14). The serial cluster of bound surfactant molecules with the intrapolymer interactions can be treated in 
the present model. A series of bound states with associating surfactants on a polyion chain is considered 
as a ‘cluster’. The cluster, whose size (number of associating ligands) is more than one, is serial b 2 states 
having bL and bR states on both ends, such as (bRb 2 ---b 2 bL). The bo state without association is also 
assumed to be the smallest cluster whose size is one. In figurethere are three clusters containing an 
isolated surfactant. Then, the fraction of the cluster Is given by [bo] -i- [bR]. The mean size of the cluster Iq 
is defined as. 


h - 


B 

[bo] + [bL] ’ 


(2.4) 


where Iq 1. By a simple idea, at a full binding state 6 = 1, the following relation is assumed 


[bL] _ [bR] _ [b2] _ [bz] 
[bo] [bo] [bL] [bR] ■ 


ft is expected that the Iq value tends to w-i-1 when 6 approaches unity. The distribution of the cluster size 
can be obtained as well. 


2.2. Model calculation 

In the MC process, the state transition of each site is decided by the probability function with the 
change of its energy using a pseudo-random number (figure |^. The free energy change of only ligand 
binding without association between bound ligands (f ^ bo) is 

Aeb = Ae“ - rin(Cf) = -A;b Tin (^TCf), (2.5) 

which is dependent on the concentration of free ligand. The free energy change with association between 
a pair of bound ligands (bo ^ bL or bR and bR or bR ^ bz) is 


Aea =-fcBrinu. (2.6) 

The acceptance of state change trials Is determined according to the standard Metropolis algorithm using 
the probability function of the energy change of the system I1SII16I . 
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For example, if the state of the site is bo, the probabilities of the MC trial for the transitioit to the states 
f, bL and bR are 1/2, 1/4, and 1/4, respectively. In the process from bo to bR, if the bound ligand which 
can associate to the objective ligand in the left-hand side within width n is found, the state transition 
is decided in the above mentioned manner. If no associable ligand is found, the process is discarded. 
Other probabilities of the process and the manner of decision are similarly given in order to maintain the 
detailed balance (see figurej^. 

The number of MC step per site is chosen within the range from 10® to 10® at a given KCf (Acb ) value 
so that the convergence of the quantities is achieved in the early stage (mostly within a few percent of 
half MC steps). After half discarded MC steps, the quantities are finally evaluated by summing up in the 
later MC steps. 

All f state is usually adopted as an initial state of sites for a series of calculations. A final state at each 
KCf value is used as an initial state of the next point. The MC trials are carried out in the order of the 
lattice site number within a single loop. The results are unaffected even if the site is randomly chosen, or 
an initial state is changed. 

In the standard calculations, the lattice has N - 1000 sites under the periodic boundary condition. 
When N - 5000, no difference of the result is observed. The generalized Fibonacci method is used as a 
random number generator (17). 


2.3. Analytical solutions 


In the classical cooperative model without various more realistic effect, we have the well-known 
equation of binding fraction by Satake and Yang (3, expressed as. 


1 

l-s 

1 

1 

^ 1 - KuCf 

2 

\/(l - s)2- 1 -4crs 

2 

^/il-KuCf]^+4KCf 


where well-used parameter s is defined as, s = KuC^ (7). When u-\, the relation means simple binding, 
as 




KCf 


(2.8) 


\ + KCf 

The binding isotherm by our previous scheme with a matrix method is also presented in the case of 
n = 2 in figure|^and figure |^(T0). 


3. Results and discussion 

3.1. Dependence on the cooperativity parameter u 

First, the fundamental properties of our model system are presented as compared with the classical 
cooperative bindings of equation |2.7| . The dependence of binding isotherm on the cooperativity param¬ 
eter u is presented at the association width n = 1 in figurej^ In this case, the association with the nearest 
neighbor is only considered. At u = 1, the MC calculation curve shows a higher affinity (lower CAC) and 
cooperativity than the simple binding curve of equation (2.8| due to the present binding scheme with four 
bound states. With an increase of the parameter u, the difference between both curves becomes smaller. 
At w = 100, the differences of Q values are within +10% at most points. 

In figure]^ the dependence on the cooperativity parameter u is shown in the case of n = 4. Associ¬ 
ations between the surfactants separated by binding sites within the width n = 4 reduce the CAC values 
to about one third of the classical model. The association with larger u value results in higher affinity 
and asymmetrical binding isotherms in the case of the ligand association across the binding sites. In the 
high u cases, most of the bound surfactants provide associations with the neighboring ones in the range 
9 < 0.5. Then, the bindings become relatively weak in high 9 range. At u = 100, the MC result shows a 
slight two-phase binding curve with a shoulder at around 9 - 0.5, corresponding to the rearrangement of 
clusters (shown below). 

The binding isotherm with a sharp rise and gradual saturation are frequently observed in the experi¬ 
ments of surfactant binding to a flexible polyion(4l|^. A factor deforming binding isotherm is exhibited 
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- 3 - 2-1012 


log(/fCf) 


Figure 3. Dependence of binding isotherm on the cooperativity parameter u at association width k = 1. 
Symbols, MC calculations, lines, curves by the Satake-Yang equation, equation u; 1,3.33 (1/m = 0.3), 
10, 33.3 (1/m = 0.03), and 100, from right to left, respectively. N = 1000 under the periodic boundary 
condition. 



-3 - 2-10 1 

logdfCf) 

Figure 4. The dependence of binding isotherm on the cooperativity parameter u at association width 
n = 4. the others are the same as in figurej^ 


from our model calculations. In the previous works, we showed that the two-phase binding isotherm is 
interpreted by the hypothesis of ligand association across a site and/or of steric hindrance. In the present 
model, we can observe the similar feature only in the restricted (n-2-8 and high u) cases. 

3.2. Dependence on the association width n 

Dependence of the binding isotherm on the association width n is shown in figure in the case of 
u - 10. With an increase of the association width n, the binding becomes stronger and the curve clearly 
shows non-symmetrical feature with a sharp rise and relatively gradual saturation. In the cases of n = 2 
or 4, cooperative binding looks steep in low 6 range, but in larger n{n^8) cases, the cooperativity of the 
isotherm becomes low. 

In a more cooperative case of u = 100, the dependence of the curves on n is complicated as shown 
in figure]^ The distances between the binding isotherms are slightly wider than those in figurej^ When 
n = 2 or 4, a strong two-phase feature appears with a shoulder on the binding isotherm. In large n(n>8) 
cases, the binding isotherm becomes monotonous with a very sharp rise. 

The shoulders of binding isotherms in figurej^and figurej^are interpreted as follows. In the cases of 
n = 2-4 at high cooperativity parameter for low binding fraction (0 < O.S), the ligands tend to entropically 
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-3 -2 -1 0 

\og{KCi) 

Figure 5. The dependence of binding isotherm on the association width n at cooperativity parameter 
u = 10. n\ 1 (x), 2 (square), 4 (triangle), 8 (cross), 16 (circle), and 32 (inverted triangle), from right to left, 
respectively. N = 1000 under the periodic boundary condition. Full line is drawn by the Satake-Yang’s 
theory. Broken line indicates the isotherm of n = 2 case calculated in our previous scheme ()Sia = /)ib = 

= 1, y= 1, (T= 1 /m = 0.1) (To). 



-4 -3 -2 -1 

logfATCf) 

Figure 6. The dependence of binding isotherm on the association width n&tu- 100. The others are the 
same as in figure)^ Broken line indicates the isotherm of « = 2 case calculated in our previous scheme 
(;6ia = /fib = /3i = 1, r = 1. o- = 1/u = 0.01). 


bind at intervals, leading to associations across unoccupied site(s). When the binding fraction becomes 
around O.S, the ligands begin to bind to the remaining unoccupied sites. However, these bindings are 
suppressed because the newly bound ligand cannot sufficiently associate with the neighboring ligands. In 
the higher binding range, recombination of the associations occurs, represented by complicated cluster 
size changes shown in figure Such a transition of the situation of the ligands makes a shoulder on 
the binding isotherms. In the cases of n = 5 - 8 at high cooperativity parameter, the binding isotherm 
frequently shows a stepwise increase of 9. In the cases of still larger n, the transition becomes obscure 
due to the extensive possibilities of association formations. 

The isotherms in our previous scheme can be calculated only in the n - 2 case using the matrix 
method (broken hne in figureand figure]^. Although the curves and CAC values resemble those of 
n = 2 in the model, their rising tendency is more gradual in high 9 range than those of the present model 
due to the incorporation of the flexibility in the present model. 

It has been observed that the long tail of the surfactant molecule enhances the binding affinity, and 
leads to low CAC (2|T8l[^. In these cases, the long surfactant may induce an increase in association 
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region (association width n in the model). This can effect the shape of the binding isotherm as well as 
the CAC value. The effect of molecular geometry of surfactant must be reconsidered from this point of 
view I11II20I . 

3.3. End effect (Dependence on the lattice length) 

The above MC results are obtained under the periodic boundary condition of AT = 1000, i.e., by con¬ 
necting both ends of the lattice. The end effect and the dependence of the polyion (lattice) length are 
examined without the periodic boundary condition. Figure shows a t)rpical effect of small N values on 
the Isotherms under non-perlodlc condition at u- 100 and n- 4, compared with N - 1000 case under 
the periodic boundary condition. 



-3 - 2.5 -2 - 1.5 -1 


log(li:Cf) 

Figure 7. The dependence of binding isotherm on the number of lattice points N at u = 100 and n = 4 
without periodic boundary condition. N; (closed symbols) 1000 (triangle), 300 (inverted triangle), 100 
(right-hand triangle, dark), 30 (right-hand triangle, light), and 10 (square), from left to right, respectively. 

Open triangle; N = 1000 case under periodic boundary condition. 

In the case of N - 1000, the end effect cannot be detected, and in the cases of N > 100 the effect is 
relatively small. In M < 100 cases, however, the shift of the binding curve is significant. The decrease in 
N makes the binding curve smooth and gradual. The end effect Is larger In larger n and/or larger u cases 
(data not shown). 

The effect of polyion size on the surfactant-polyion system has been investigated by several research 
groups 12111251 . With a decrease of the polymer chain length, an increase of CAC and a decrease of coop- 
erativity have been observed |2Tl|22l. The reduction of the binding with a decrease of the polyion length 
can be interpreted by this calculation. 

3.4. Cluster size 

In figure the changes of mean cluster size with the surfactant binding are shown with the same 
parameters as in figure The onset of self-assembly Is coincident with the CAC depending on n value. 
In cases of n-2, 4, and 8, the curves show up and down in a complex manner. The cluster size distribu¬ 
tions are monotonous due to the one-dimensional lattice system (data not shown). In cases of n = 2 - 8, 
ligands separated on the lattice induce aggregations in low 9 range with collapse of the polylon chain. 
The rearrangements of the association between ligands occur with an increase of the ligand binding in 
the middle 6 range. In cases of larger n, the rearrangement is more gradual and spread to wider 6 range. 
This unexpected behavior may be important to understand the binding of the surfactant to the flexible 
polylon. 

In the model calculation, large u value leads to high /c- h has been shown that a large cluster size 
(aggregation number) corresponds to a high cooperativity (26). The dependence of the cooperativity on 
the charge density of the polyion may be interpreted in the present idea. 
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logiKCt) 


Figure 8. The dependence of variation in mean cluster size Ic on the association width n at cooperativity 
parameter u = 100. Symbols are the same as in figurej^ 


In the present model, one ligand can simply associate with two ligands on both sides of the polyion 
chain. More complicated models showing phase transition, should be suitable for the coil-globule transi¬ 
tion by introducing additional hypotheses I27II28I . 

3.5. Comparison with the experiment 

Shimizu reported the bindings of cationic surfactants by several flexible poly(carboxylic acid)s using 
a potentiometric technique (29|. Most of the binding isotherms show a sharp onset followed by gradual 
saturation. The gradual slope in a range of higher 6 means that highly ligand association suppresses a 
new binding, although the interpretation due to the electrostatic factor is also likely. 

The tentative comparison of the calculated binding curve with an experimental one is executable, 
although the precise data-fitting in the whole range is difficult. Figure shows the binding isotherm of 
the dodecylpyridinium ion (DP+) by poly acrylic acid (PAA) for a - 1.0 at 0.01 m NaCl (29) with the corre¬ 
spondent calculated curve. In the low surfactant concentration range, reasonable agreement is achieved 
with the association width n-3 and estimated cooperativity parameter a -Hu- 0.003. The high coop- 



- 4.5 -4 - 3.5 -3 - 2.5 


log(Cf) 

Figure 9. Comparison of the experimental isotherm of dodecylpyridinium ion to poly acrylic acid for 
a = 1.0 (large squares) by Shimizu (29) with the correspondent calculated values (smaU filled squares). 
Experimental conditions; T - 30 °C, = 0.01 mol-kg“^ H 2 O, and mp = 0.915 mN. Parameters of 

the model calculation; N = 1000, n - 3, a = llu = 0.003, and logff = 0.66 = -3.8 kj/mol, Aca = 

-14.6kJ/mol). 
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erativity and the association within the third neighboring ligands lead to this binding profile. The slight 
shoulder of the experimental isotherm in the higher range may be due to additional weak associations of 
the ligands, as observed in other isotherms (29). 


4. Conclusions 

The objective of the study is to analyze the essential characteristics of the system using a simple 
model. It seems more suitable rather than molecular simulation studies. The present model is examined 
to investigate the binding mechanism of ionic surfactants to a flexible polyion in dilute solution by the MC 
calculation. Introducing the non-cooperative binding and the ligand association separated by multiple 
lattice points, we can show their effects as new factors that determine the binding isotherm. Charges 
on the polyion should play an appreciable role in increasing the local concentration of the surfactant 
molecules. 

In this study, a simple model is examined by the MC calculation to investigate the binding mechanism 
of ionic surfactants to a flexible polyion in dilute solution. Charges on the polyion play a role in increas¬ 
ing the local concentration of the surfactant molecules. Introducing the ligand association separated by 
lattice points within the association width n in the stochastic processes, we can show its effect as a new 
factor that determines the binding isotherm. 

In the case at high uoi n- 1, the results of the present model converge to those of a classical coop¬ 
erative model. With making the width n double, the binding affinity looks approximately twofold, i.e., 
reducing CAC by half. In the range of n = 2 - 8, two-phase or stepwise nature of the isotherm is observed 
due to the complicated rearrangement of the surfactant cluster on the polylon chain, in particular, in the 
high u cases. The short polyion chain also lowers the affinity and their cooperativity In the model due to 
the end effect. 

The model of cluster formation is so simple that the results are rather limited at the present stage. 
However, our approach suggests the essential points of the effect of polyion flexibility on the binding 
isotherm reflecting the freedom of ligand binding and of their association. 


Acknowledgements 

The calculations in this study were carried out using the workstation in the information technology 
center of the Hamamatsu University School of Medicine. 


References 

1. Polymer-Surfactant Systems. Surfactant science series, Voi. 77, Kwak J.C.T. (Ed.), Marcei Dekker, New York, 1998. 

2. Kogej K., Adv. CoUoid Interface Sci., 2010,158, 68; doi 10.1016/].cis.2009.04.003 

3. Chiappisi L., Hoffmann L, Gradzielski M., Soft Matter, 2013, 9, 3896; doi: 10.1039/C3SM27698H 

4. Shimizu T., In: Polymer Materials Encyclopedia, Voi. 8, Salamone J.C. (Ed.), CRC Press, Boca Raton, 1996, 5783- 
5791. 

5. Mata J., Patel J., Jain N., Ghosh G., Bahadur P, J. CoUoid Interface Sci., 2006,297,797; doi 10.1016/].icis.2005.11.022 

6. Schwarz G., Eur. J. Biochem., 1970,12, 442; doi 10.1111/i.l432-1033.1970.tb00871.x 

7. Satake L, Yang J.T., Biopoiymers, 1976,15, 2263; doi 10.1002/bip.1976.360151113 

8. Nishio T., Shimizu T., Kwak J.C.T., Minakata A., Biophys. Chem., 2003,104, 501; 
doi: 10.1016/50301-4622(03)00039-5 

9. Yan P, Jin C., Wang C., Ye J., Xiao J.-X., J. CoUoid Interface Sci., 2005, 282,188; doi 10.1016/j.jcis.2004.08.164 

10. Nishio T., Shimizu T., Biophys. Chem., 2005,117,19; doi 10.1016/j.bpc.2005.03.011 

11. Jain N., Trabelsi S., Guillot S., McLoughlin D., Langevin D., LeteUier P, Turmine M., Langmuir, 2004, 20, 8496; 
doi:10.1021Aa0489918 

12. Lapitsky Y, ParUth M., Kaler E.W., J. Phys. Chem. B, 2007, 111, 8379; doi: 10.1021/jp0678958 

13. Dias R.S., Magno L.M., Valente A.J.M., Das D., Das P.K., Maiti S., Miguel M.G., Lindman B., J. Phys. Chem. B, 2008, 
112, 14446; doi 10.1021/ip8027935 


43302-10 




Binding of ionic surfactants to fiexibie poiyions 


14. Nizri G., Lagerge S., Kamyshny A., Major D.T., Magdassi S., J. Colloid Interface Sci., 2008, 320, 74; 
doi: 10.1016/j.jcis.2008.01.016 

15. Metropolis N., Rosenbluth A.W., Rosenbluth M.N., Teller A.H., Teller E., J. Chem. Phys., 1953, 21,1087; 
doi: 10.1063/1.1699114 

16. Nishio T., Biophys. Chem., 1991,40,19; doll0.1016/0301-4622(91)85026-M 

17. Brent R.P., In: Proceedings of Conference “Fifth Australian Supercomputer Conference” (Melbourne, 1992), Mel¬ 
bourne, 1992, 95-104. 

18. Ishiguro M., KoopalL.K., CoUoids Surf A, 2009, 347, 69; doi:10.1016/j.colsurfa.2008.12.008 

19. Liu J., Nakama M., Takisawa N., Shirahama K., Colloids Surf A, 1999,150, 275; 
doi:10.1016/S0927-7757(98)00819-X 

20. Malovikova A., Hayakawa K., Kwak J.C.T., J. Phys. Chem., 1984, 88,1930; doi 10.1021/jl50654a002 

21. Liu J., Takisawa N., Shirahama K., Abe H., Sakamoto K., J. Phys. Chem. B, 1997,101, 7520; doi 10.1021/jp9711981 

22. Liu J., Shirahama K., Miyajima T, Kwak J.C.T., Colloid Polym. Scf, 1998, 276,40; doi: 10.1007/s003960050206 

23. Shimizu T, J. Phys. Chem. B, 2003,107, 8228; doi 10.1021/jp022541b 

24. Meszaros R., Varga I., Gilanyi T., J. Phys. Chem. B, 2005,109,13538; doi 10.1021/jp051272x 

25. D’Errico G., CiccareUi D., Ortona O., Paduano L., Sartorio R., J. CoUoid Interface Sci., 2007, 314, 242; 
doi: 10.1016/j.jcis.2007.05.035 

26. Hansson P., Almgren M., J. Phys. Chem., 1996,100, 9038; doi 10.1021/jp953637r 

27. Mel’nikov S.M., Sergeyev V.G., Yoshikawa K., J. Am. Chem. Soc., 1995,117, 9951; doi 10.1021/ia00145a003 

28. Kuhn P.S., Diehl A., Phys. Rev. E, 2007, 76, 041807; doi 10.1103/PhysRevE.76.041807 

29. Shimizu T, Colloids Surf A, 1995, 94,115; doi 10.1016/0927-7757(94)02986-5 


MoAe/ibHe AOC/ii/pKeHHJi KO/ieKTi/iBHoro SB'nsyBaHHB Iohhi/ix 
cyp(|)aKTaHTiB 3 npoTi/i;ie>KHO 3apjiA>KeHi/iMi/i 
THyMKi/iMi/i no/iiioHaMi/i 

T. HiiiiicP, T. llliMis^, HI. MoiiiiflcPI A. MiHaKaTi/P 

^ Biflflm o6'eflHaHMX Hayx npo arofli/iHy (((lisuKa), Meflunna ujKO/ia yHiBepcMrery m. XaMapy, 

XaMapy 431-3192, Bnonin 

^ Biflflm iHxeHepn e.fieKTpoHHnx ra iHtjjopMapiilHi/ix cncreM, (jjaxy/ibrer npi/ipoflHUHi/ix nayx i rexHoaorii, 
yniBepcMTeT m. Xipocaxi, Xipocaxi 036-8561, Bnohin 


/(/in BMBHeHHB BsacMOAM MO/iexy/i iOHHOro cyptfiaKTaHTa 3 rnyuKnM no/iiiOHHi/iM /lanpio/KKOM y posBefleno- 
My po34MHi 3anponoHOBaHO noBy CTarncTi/iHHy MOfle/ib A/in KO/ieKTi/iBHOro 3B'n3yBaHHn MOHOMepHnx /liranfl 
3 /liHifiHOto rpaTKOto. flpi/inycKacTbcn, u^o e/ieKipocTari/iuHe 3B'n3yBaHHn /liranfli/i 3 By3/iOM na no/iiiOHi ra ri- 
flpoc()o6Hi 3B'n3Ki/i Mix cyciflHiMi 3B'n3aHi/iMn /liraHflaMi/i e croxacn/iHHi/iMM npopecaMi/i. /l/in rHytKoro no/iiiona 
BBOflUTbcn acopiapin /liranfl, BiflOKpeM/ienux flexi/ibKOMa rparKOBMMi/i TOuxaMn b Mexax Bi/iaHaneHOi uji/ipnHi/i. 
/(/in Toro, i±(o6 floc/iifluri/i noBeflinxy 3B'n3yBaHHn, afliOcHOtOTbcn o6HMC/ieHHn MeroflOM Monre Kap/io. rinoreaa 
acopiapii/liraHflu ra iT lunpi/iHn na /lanptoxKy e Ba>K/ii/iBOto fl/in Bi/i3HaHeHHn KpnmuHOi KonpeHTpapii arperapii 
ra i30TepMi/i 3B'n3yBaHHn. Peay/ibra™ e npnilHnTHMMM fl/in iHrepopernpiT flexi/ibKOx excnepi/iMeHTiB no 3B'n3y- 
BaHHto cyptjjaKTaHT-rHynKi/ifi no/iiiOH. flpeflcraB/ieHO ra oOroBopeno aacrocyBaHnn Merofly. 

KniOMOBi c/iOBa: ssae/woyjw /ohh/im cypcjiaKTaHT-my^Kua no/iiiOH, rpaiKa niHiUHOro noniiona, KoneKmane 
3B'ft3yBaHHf! /liraHfl, boTepMi/t 3B'B3yBaHHB, LUMpi/iHa acopiapif, MOflejibne oSm^cnenHa 
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